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ABSTRACT 

We present the first far ultraviolet (FUV) spectra of the four known Balmer-dominated supernova 
remnants (SNRs) in the Large Magellanic Cloud, acquired with the Far Ultraviolet Spectroscopic 
Explorer {FUSE). The remnants DEM L 71 (SNR 0505-67.9), SNR 0509-67.5, SNR0519-69.0 and 
SNR 0548—70.4 arc all in the non-radiative stages of evolution and exhibit expansion speeds ranging 
from ~ 500 kms _1 to ~ 5000 kms -1 . We have detected broad emission lines of Ly (3, Ly 7, C III 
and O VI in DEM L 71 (Vfwhm ~ 1000 kms -1 ) and have detected broad Ly (3 and O VI emission in 
SNR 0519-69.0 {Vfwhm -3000 kms" 1 ). In addition, broad Ly f3 emission (Vfwhm ~3700 kms" 1 ) 
has been observed in SNR 0509—67.5, the first detection of broad line emission from this SNR. No emis- 
sion was detected in our FUSE spectrum of SNR 0548—70.4, allowing us to place only upper limits 
on the FUV line fluxes. The spectra of these SNRs are unaffected by postshock cooling, and provide 
valuable probes of collisionless heating efficiency in high Mach number shocks. We have used the F(Ly 
0) /Fo v/(1032) flux ratio and relative widths of the broad Ly (3 and O VI lines to estimate the degree of 
electron-proton and proton-oxygen ion equilibration in DEM L 71, SNR 0509—67.5 and SNR 0519—69.0. 
Although our equilibration estimates are subject to considerable uncertainty due to the faintness of the 
FUV lines and contributions from bulk Doppler broadening, our results are consistent with a declining 
efficiency of electron-proton and proton-oxygen ion equilibration with increasing shock speed. From 
our shock velocity estimates we obtain ages of 295-585 years for SNR 0509—67.5 and 520-900 years for 
SNR 0519— 69.0, respectively, in good agreement with the ages obtained from SN light echo studies. 

Subject headings: ISM: supernova remnants: individual (0505-67.9 (DEM L 71), SNR 0509-67.5, 
SNR 0519-69.0, SNR 0548-70.4)- ISM: kinematics and dynamics, shock waves 



1. INTRODUCTION 

The proximity of the Large Magellanic Cloud (LMC) 
(50 kpc; Feast 1991, 1999) and the relatively small extinc- 
tion along the line of sight to this galaxy (E(B-V) ~ 0.1) 
make it an excellent laboratory for study in the X-ray and 
ultraviolet. In particular, the well calibrated distance of 
the LMC allows us to study both the kinematics of indi- 
vidual supernova remnants (SNRs) (Blair et al. 2006) 
and properties of the global SNR distribution within this 
galaxy. 

The first systematic X-ray survey of the LMC was per- 
formed with the Imaging Proportional Counter of the Ein- 
stein observatory (Long, Hclfand & Grabelsky 1981) in the 
0.15-4.5 keV range. The observations revealed that nearly 
half of the 97 detected sources exhibited an extended 
morphology characteristic of SNRs, making these objects 
the largest class of X-ray sources revealed in the survey. 
Four of the SNRs identified with Einstein were shown by 
Tuohy et al. (1982) to exhibit optical spectra consist- 
ing almost entirely of faint hydrogen Balmer line emission 
(i.e., they are 'Balmer-dominated'). The four remnants — 
0505-67.9 (DEM L 71 from the Davies, Elliott & Meaburn 
(1971) catalogue), SNR 0509-67.5, SNR 0519-69.0 and 
SNR 0548—70.4 vary widely in age, ranging from a mere 



few hundred years (SNR 0509-67.5; Rest et al. 2005) to 
-10,000 years (SNR 0548-70.4, Smith et al. 1991). 

Balmer-dominated line emission is produced when the 
expanding SNR encounters low density (—0.1-1 cm -3 ), 
partially neutral interstellar ISM (Chevalier & Raymond 
1978; Bychkov & Lebedev 1979). Due to the high tem- 
perature and low density behind the blast wave, the cool- 
ing time scale of the gas greatly exceeds the dynamical 
time scale of these SNRs (they are non-radiative). The 
hydrogen overrun by the shock is collisionally excited in a 
thin (<10 15 cm) ionization zone, producing optical spectra 
dominated by the Balmer lines of hydrogen (Chevalier & 
Raymond 1978; Chevalier, Kirshner & Raymond 1980). 
In the ultraviolet, the spectra are 'Lyman-dominated'. 
Balmer-dominated shocks are mostly found in young rem- 
nants (e.g., SNRs in the Sedov- Taylor stage or younger) of 
Type la explosions. In these objects the survival of signif- 
icantly neutral, low density ISM gas is likely aided by the 
lack of strong ionizing radiation and stellar wind from the 
progenitor star (a white dwarf). X-ray spectroscopy of the 
reverse-shocked ejecta in the four LMC Balmer-dominated 
remnants has revealed strong Fe L shell emission (Hughes 
et al. 2003; Hendrick, Borkowski & Reynolds 2003; War- 
ren & Hughes 2004), confirming the Type la origin of 
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these SNRs. These overall characteristics place the four 
Balmer-dominated LMC remnants into the same class as 
the well-known Galactic Type la SNRs Tycho and SN 1006 
(Chevalier, Kirshner & Raymond 1980; Kirshner, Winkler 
& Chevalier 1987; Smith et al. 1991; Ghavamian ct al. 
2001), making them the first examples of such remnants 
identified outside our galaxy. 

The hydrogen line emission in Balmer-dominated shocks 
consists of two components: a narrow component pro- 
duced by collisional excitation of overrun ambient H I, 
and a broad component generated when the hot postshock 
H I produced by charge exchange is collisionally excited. 
Both broad and narrow component Balmer lines can often 
be detected in the optical, but interstellar absorption in 
the ultraviolet completely removes the narrow component 
Lyman line emission. However, if the broad component 
Lyman lines are both bright and wide, the profile wings 
can be detected. 

The diagnostic utility of the broad H I lines is twofold. 
First, the width of a broad H I line gives the postshock pro- 
ton temperature, thus yielding an estimate of the shock 
speed when combined with the Rankinc-Hugoniot jump 
conditions (Chevalier, Kirshner & Raymond 1980; Smith 
et al. 1991; Ghavamian et al. 2001). Second, comparing 
the width and flux of the Lyman lines with those of ultra- 
violet resonance lines such as O VI AA1032, 1038 and C III 
A977 yields an estimate of the ion-ion temperature equi- 
libration at the shock front, well before the equipartition 
of energy between the different ionic species by Coulomb 
collisions (Hester, Raymond & Blair 1994; Raymond, Blair 
& Long 1995). This is key information for probing the na- 
ture of collisionless (plasma) heating processes at the shock 
front, one of the central problems in understanding high 
Mach number shocks in many astrophysical settings. 

The availability of broad hydrogen, carbon and oxygen 
emission lines, along with high quality X-ray spectra from 
Chandra provides us with an invaluable opportunity to si- 
multaneously study both the electron-proton and ion-ion 
equilibration in fast collisionless shocks and the evolution 
of remnants from Type la SNe. In this paper we present 
FUV spectra of the four Balmer-dominated LMC SNRs. 
We utilize numerical shock models to interpret the spectra, 
constrain the kinematics of the four SNRs and investigate 
collisionless heating processes in fast astrophysical shock 
waves. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. THE RAW DATA 

We observed the four Balmer-dominated LMC SNRs 
with the Far Ultraviolet Spectroscopic Explorer (FUSE) 
in 2001 September 23 and 24 (see Moos et al. 2000 
and Sahnow ct al. 2000 for a detailed description of the 
spectrograph). The data were acquired as part of a Cy- 
cle 3 Guaranteed Time Observation program (P214; W. 
P. Blair, PI), with each remnant observed through the low 
resolution (LWRS) square aperture of FUSE (30" x 30"). 
In each case the LWRS aperture was centered on the ge- 
ometric center of each SNR as measured in the optical by 
Tuohy et al. (1982). The roll angles of these observa- 
tions were set by the nominal spacecraft roll at the time 
of data acquisition. The total integration time and the 

6 http: / /fuse. pha.jhu.edu/support/guidc/guide. html 



time obtained during orbital night is listed for each object 
in Table 1. The Balmer-dominated LMC SNRs vary in 
size from approximately 30" across (SNR 0509—67.5) to 
l.'l across (SNR 0548—70.4), resulting in significant varia- 
tions in the fraction of the shell covered in each observation 
(see Figure 1). While the LWRS aperture covered nearly 
the entire shell of SNR 0509-67.5 and SNR 0519-69.0, it 
sampled primarily the approaching and receding faces of 
the shell in DEM L 71 and SNR 0548-70.4. As we show in 
the following sections, the fraction of each SNR covered by 
the FUSE aperture can significantly affect the shapes and 
centroids of the FUV lines in the final extracted spectra. 

We acquired a separate FUSE observation of DEM L71 
during Cycle 3 (Guest Observer program C072; PI: P. 
Ghavamian). In this observation we placed the MDRS 
aperture (4" x20") parallel to the southern edge of the 
blast wave at a position angle of 113° (see Figure 1). The 
aperture orientations were nearly constant during these 
observations. A total integration time of 121.6 ks was 
accumulated during five separate pointings on 2003 Jan- 
uary 4-7 (C0720202-C0720205) and 2003 December 29 
(C0720206). However, the target acquisition in observa- 
tion C0720205 failed, leaving 100 ks of useful data. 

Our data reduction was performed using version 3.1.3 of 
the CalFUSE pipeline. The expected FUV line emission 
in the four LMC remnants is expected to be faint (<^10~ 13 
ergs cm~ 2 s _1 ), with line widths (j>1000 kms -1 ) greatly 
exceeding the resolution of the filled LWRS aperture (ap- 
proximately 106 kms -1 ). These properties make the de- 
tection of the FUV lines quite sensitive to the quality of 
the background subtraction. Since the CalFUSE pipeline 
derives an empirical background for each spectrum us- 
ing non-illuminated portions of the detector, we can im- 
prove the signal-to-noise of the measured background by 
co-adding the individual raw exposures for each of the four 
FUSE detectors (1A,1B,2A,2B), then extracting 1-D spec- 
tra from the combined raw frames. We used the IDL rou- 
tine ttag_combine to combine the raw exposures (Figure 2), 
then used the CalFUSE pipeline to extract the 1-D spec- 
tra. Wc retained only the orbital night portions from each 
dataset (Table 1) to minimize the contribution of airglow 
emission lines and scattered airglow continuum to the ob- 
ject spectra. 

We combined the raw exposures from the MDRS data 
in the same manner described above for the LWRS data. 
Since the C0720206 data were acquired nearly a year after 
the other observations, the FUV spectrum of this obser- 
vation is significantly shifted in detector pixels relative to 
the earlier observations. Therefore, we first co-added the 
C0720202-C0720204 datascts into one raw image, then co- 
added C0720206 into a separate raw image. We extracted 
spectra from each image, cross-correlated each one and fi- 
nally co-added them to obtain the 1-D MDRS spectrum 
of the DEM L71 blast wave. 

2.2. SPECTRAL EXTRACTION 

The FUSE spectra are recorded on four channels: LiF 
1, LiF2, SiC 1 and SiC 2. Each channel is divided into 
A and B segments, with the two segments covering the 
approximate wavelength range 905 A - 1100 A for each of 
the LiF 1 and LiF 2 channels, and the wavelength range 
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987 A- 1180 A for the SiC 1 and SiC 2 channels (see 
the FUSE Observer's Guide 6 ). In each of our LWRS ob- 
servations, the LiF 1A and SiC 2A (917 A- 1006 A) seg- 
ments offer the highest effective area at the positions of 
the expected FUV emission lines. Therefore, we searched 
the spectra from these segments for SNR emission first. 
Although the extracted FUSE spectra are meant to be 
summed into one high signal-to- noise data product, inter- 
channel thermal drifts during each exposure can result in 
misalignments as large as 10" between LiF 1 and the other 
channels. While the LiF 1 channel was held fixed on the 
SNR during an observation, the others sampled emission 
from slightly different portions of the SNR. The result- 
ing spectra can differ from channel to channel. Therefore, 
for each observed SNR we combined spectra from multi- 
ple channels only when they appeared consistent with one 
another to within the flux uncertainties. 

In the case of DEM L 71, prominent, broad emission 
lines of Ly (3, Ly 7, O VI and C III are detected in the 
LWRS observations. The spectral features are observed 
in all segments of the FUSE detector save for the LiF IB 
and LiF 2A segments, whose spectral range excludes these 
emission lines. Although inter-channel drift is expected to 
have occurred during the 28 ks observation, the spectra are 
all consistent with one another to within the errors. We 
generated the spectra in Figure 3 by cross-correlating and 
combining data from SiC 2A with SiC IB (905 A - 1005 
A) and the data from LiF 1A with those of LiF 2A and 
SiC 1A (988 A - 1090 A). Although broad Ly (3 and O VI 
emission were seen in the SiC 2B segment, we did not in- 
clude this data in our co-addition because the background 
was unevenly over-subtracted for this segment by the Cal- 
FUSE pipeline. 

Unlike DEM L 71, SNR 0509-67.5 and SNR 0519-69.0 
fit almost completely within the LWRS aperture (Fig- 
ure 1), increasing the likelihood that inter-channel drifts 
would adversely affect the detection of FUV emission from 
these objects. Since SNR 0509—67.5 exhibits the faintest 
FUV emission of all the Balmer-dominated SNRs (Fig- 
ure 4), we restricted our analysis of this object to data 
from the LiF 1A segment. 

In the case of SNR 0519-69.0, on the other hand, FUV 
line emission was detected in all segments shortward of 
1100 A. One factor complicating our interpretation is the 
presence of a very faint continuum under the Ly (3 and 
O VI features (clearly seen in Figure 2). The continuum is 
absent on the short wavelength side of the broad Ly (3 line, 
but rises to a nearly constant level beneath the broad lines 
and continues to the long wavelength end of the spectrum. 
The continuum is absent in the other channel spectra near 
1030 A, and appears to be stellar contamination. After 
comparing the emission in all channels, we cross-correlated 
and combined the LiF 1A, LiF 2B, SiC 1A and SiC 2B data 
to obtain spectra in the 940 A— 1050 A range. Although 
broad Ly 7 was also detected, it was present solely in the 
SiC 2A spectrum of SNR 0519-69.0 (the channel with the 
highest effective area in this wavelength region). 

In our MDRS observation of DEM L 71, FUV emission 
was detected only in the spectrum of the LiF 1A segment. 
The lack of detection in the other segments is likely due 

7 IRAF is distributed by the National Optical Astronomy Observatories, 
with the National Science Foundation. 



to a combination of their lower effective areas and the the 
strong impact of inter-channel thermal drifts. The latter 
effect is especially important here, where the smaller size 
of the MDRS aperture makes it easier for non-primary 
channels to drift off the targeted filament. 

One of the more surprising results of our observations is 
the lack of any discernible FUV emission from the largest 
of the Balmer-dominated SNRs, SNR 0548-70.4 (see Fig- 
ure 6). No emission is detected in either the raw combined 
data, the individual channel spectra or the co-added spec- 
tra from all channels. In the following section we will 
describe a possible explanation for this non-detection, and 
discuss its implications. 

3. FAR ULTRAVIOLET SPECTRA 

Of the observed Balmer-dominated remnants, two ob- 
jects - DEM L 71 and SNR 0519-69.0 - exhibit read- 
ily identifiable, broad FUV line emission consistent with 
Balmer-dominated shocks: Ly (3 (A1025.7), O VI AA1032, 
1038, Ly 7 (A972.5) and C III A977. In a third remnant, 
SNR 0509-67.5, only the broad Ly (3 line is detected with 
reasonable certainty. However, this is the first detection 
of broad hydrogen line emission from this remnant. Sur- 
prisingly, no discernible FUV emission is seen from the 
fourth remnant, SNR 0548—70.4, despite the known de- 
tection of broad Balmcr line emission in this remnant by 
Smith et al. (1991). Here, we describe the FUSE spec- 
tra of each remnant in detail. All spectral fits have been 
performed using the SPECFIT task (Kriss 1994) from the 
STSDAS contributed package in IRAF 7 unless noted oth- 
erwise. 

3.1. DEM L 71 



3.1.1. LWRS Spectra 

The combined LWRS spectrum of DEM L 71 (Figure 3) 
provides an especially impressive example of the FUV 
emission from a Balmer-dominated SNR. The emission 
lines are very broad (^1000 kms -1 ), bright and clearly 
impacted by Galactic and LMC halo absorption. The pres- 
ence of strong, broad Ly (3 indicates that a significant frac- 
tion of the gas entering the blast wave is neutral, while the 
relative faintness of the C III line relative to O VI indicates 
that postshock cooling is negligible (Raymond 2001). The 
wavelengths of prominent absorption features from both 
the Milky Way and LMC are marked in Figure 3. The 
broad Ly f3, O VI A1032 and O VI A1038 lines are moder- 
ately blended, as are the Ly 7 and C III A977 lines. 

The LWRS aperture samples emission near the pro- 
jected center of DEM L 71 (Figure 2). Therefore, each 
emission line should appear doubled by the bulk motion 
of the postshock gas (± | v s ) on near and far sides of 
the SNR. The shock velocities in DEM L 71 have been 
estimated to lie between 600 kms -1 and 1100 kms -1 
(Smith et al. 1991; Ghavamian ct al. 2003; Rakowski, 
Ghavamian & Hughes 2003), so the emission lines in 
DEM L 71 should appear blue/red-shifted by 400-800 
kms -1 from the rest velocity of the LMC (which we take 
to be +275 kms -1 ). However, the presence of strong in- 
terstellar absorption near the center of each emission line 
in Figure 3 makes it difficult to distinguish whether each 

which is operated by the AURA, Inc. under cooperative agreement 
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emission line is doubled. As we show in Section 4, the line 
profile shapes observed in DEM L 71 are in fact consistent 
with the presence of two Doppler-shifted components. We 
have devised a spherical shell model of the blast wave to 
address the significant interplay between bulk and thermal 
Dopplcr motions in this spectrum. We present this model 
in the next section. 

Although a spherical shell model is required for a kine- 
matic interpretation of the spectrum in Figure 3, we can 
still estimate the observed flux in each of the broad emis- 
sion features by fitting each with a single Gaussian line. 
Prior to fitting each spectrum we estimated the baseline 
level by fitting the data on the blue and red ends of the 
spectrum (devoid of emission features) with a straight line 
of zero slope. After obtaining the best fit value to the 
baseline, we held this quantity fixed during our subsequent 
fits to the emission lines. We chose a 9-component fit for 
the DEM L 71 spectrum, fitting the emission components 
with a flat baseline and three Gaussian emission lines for 
broad Ly (3 and O VI AA1032, 1038. We included three 
absorption components for broad Ly (3 (Galactic and LMC 
halo H I absorption, Galactic O I absorption at 1027 A) 
while excluding the wavelength range covering the Ly (3 
airglow line. We also included one O VI LMC absorption 
component for each of the broad O VI lines. Although 
there is possible evidence for more absorption components 
along the line of sight, we found that adding more of these 
components did not substantially improve the statistical 
quality of the overall fit. 

We improved the constraining power of our fits by ty- 
ing together the widths and velocity centroids of the O VI 
AA1032, 1038 emission lines. We applied the same con- 
straint to the O VI absorption features, with the additional 
requirement T1032 : Ti038 = 2:1. In fitting the emission near 
972 A we used a 5 component model consisting of a linear 
baseline, one Gaussian each for Ly 7 and C III A977, and 
one Galactic and one LMC H I absorption component. We 
held the velocity width and line centroid of the Ly 7 line 
equal to the values measured from our Ly (3 profile fit. We 
included the H I absorption components solely as an aid for 
better fitting of the broad Ly 7 emission line. The fitted 
H I absorption lines do not actually provide reliable infor- 
mation on the H I column density because the data quality 
near the cores of the absorption profiles are heavily com- 
promised by the bright airglow Ly (3 line. Furthermore, 
we used a Gaussian profile shape to describe what are ac- 
tually saturated H I absorption lines, further invalidating 
a physical interpretation of the fitted parameters for these 
components. 

The resulting emission line fit for DEM L 71 is marked 
in Figure 3 and the predicted emission line parameters are 
listed in Table 2. The Ly (3 line width (1135±30 kms" 1 ) 
is significantly larger than that of the O VI A 1032 line 
(740±45 kms -1 ), already suggesting that the protons and 
oxygen ions have undergone significant temperature equi- 
libration at the shock front. However, we defer the task 
of estimating the specific shock velocities and comparing 
the results to existing optical measurements to Section 4. 
We note the slight velocity shift for the fitted emission line 
centroids: taking a systemic velocity of +275 kms -1 for 
the LMC, the Ly j3 and O VI A1032 lines are redshifted by 
100-150 kms" 1 (Table 2). This feature suggests that the 



emission from the far side of the DEM L 71 shell is slightly 
brighter than the near side. In Section 4 we quantify this 
result with more detailed, models of the blast wave shell. 

The prominent absorption features present in the O VI 
AA1032, 1038 lines are centered at 250±50 kms -1 , indi- 
cating that they arise in the LMC halo. The Galactic 
absorption component is considerably weaker. Our model 
fits indicate T1032 (Gal) = 1.7±0.6, with a =50±22kms- 1 
after deconvolving the measured FWHM with the filled slit 
spectral response of 106 km s _1 . The implied O VI column 
density is approximately (3.5±2.2)x 10 14 cm~ 2 . Both the 
line width and column density lie at the upper end of, but 
are fully consistent with, values measured by Howk et al. 
(2002) in their analysis of FUSE spectra of 12 early-type 
stars in the LMC. 

Aside from the O VI doublet, C III A977 in DEM L 71 is 
the only other heavy ion line detected in our LWRS spec- 
tra of the Balmer-dominated LMC remnants. Our single 
Gaussian fit (Table 2) indicates that the C III line is less 
than half as broad as the Ly (3 line, indicating that signif- 
icant temperature equilibration has taken place between 
protons and carbon ions at the shock front. 

3.1.2. MDRS Spectrum 

The night-only MDRS spectrum of DEM L 71 (Fig- 
ure 7) reveals prominent, broad Ly (3 emission moderately 
blended with faint O VI lines. As with the face-on LWRS 
spectrum, we fit the three lines together while including a 
fixed baseline level and both Galactic and LMC H I ab- 
sorption. The O VI lines arc faint enough to make the 
Galactic and LMC halo absorption features undiscernible, 
so we did not include absorption in our O VI profile fits. 
The inclusion of absorption components had minimal im- 
pact on the quality of the fit. The final results, quoted in 
Table 2, indicate a broad Ly (3 width of 1365±75 kms -1 
FWHM. The O VI line is significantly narrower, 935±125 
km s _1 , an indication that partial ion- ion equilibration has 
taken place at the shock transition. 

3.2. SNR0509-67.5 

In their optical spectroscopy of this SNR Tuohy ct al. 
(1982) and Smith et al. (1991) were unable to detect 
broad Balmer line emission and concluded that the line 
was too faint and broad to be detected. Our LWRS 
spectrum of SNR 0509—67.5 reveals a faint, strongly 
blucshifted broad Ly (3 line (Figure 4), the first broad emis- 
sion line detected from this SNR to date. The detected line 
emission is indeed faint, amounting to only 20% and 35% 
of the broad Ly (3 flux in DEM L 71 and SNR 0519-69.0, 
respectively. We fit the broad Ly (3 profile with a sin- 
gle Gaussian and ignored the regions containing the Ly (3 
airglow and H I absorption features. The result of this 
fit (Table 2) is a broad Ly (3 width of 3710±400 kms" 1 , 
where the uncertainty includes both the random statistical 
error (±325 kms -1 ) and the baseline uncertainty (±240 
kms -1 ; this quantity is the variation in line width when 
the background level is fixed to the upper or lower limits 
allowed by the initial baseline fit). This is the largest H I 
line width measured for a broad component in any of the 
Balmer-dominated SNRs thus far. 

Aside from being very broad, the Ly (3 line is blueshiftcd 
to — 263±125 kms -1 , indicating a large velocity shift 
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(-540 kms *) relative to the LMC frame. The veloc- 
ity shift of the broad line center is proportional to the 
projected component of the bulk velocity along the line 
of sight. The blueshift of the line indicates that the Ly 
(3 emission arises on the near side of SNR 0509—67.5. As 
seen in Figure 1, there is a prominent patch of enhanced 
Ha emission stretching inward from the southwest corner 
of the shell. It is likely that this patch lies on the near side 
of SNR 0509-67.5 and that FUSE has detected the FUV 
emission from this region. As we will show in Section 4, 
the size of the Ly (3 broad component shift is fully consis- 
tent with the position of the patchy enhancement relative 
to the geometric center of SNR 0509—67.5. Note that the 
region of enhanced emission is isolated close to the edge of 
the shell, so the range of viewing angles to this section of 
the blast wave is likely small (i.e., the emitting region can 
be approximated as a plane-parallel shock tilted slightly 
toward the observer). 

3.3. SNR 0519-69.0 

This Balmer-dominated SNR is similar in size to 
SNR 0509—67.5 (and likely of comparable age as mea- 
sured from light echoes from the SN explosion, Rest et al. 
2005). The FUV spectrum (Figure 5) exhibits broad Ly j3 
emission line at 1026 A, blended with broad O VI AA1032, 
1038 line emission. A broad Ly 7 line can clearly be dis- 
cerned near 972 A although the C III line falls below the 
detection threshold. 

We can draw a direct comparison between the broad 
Ly (3 widths measured in SNR 0519—69.0 and the Ha 
widths measured by Tuohy et al. (1982) and Smith 
et al. (1991). After deconvolving with the filled slit reso- 
lution (106 kms -1 ) we find good agreement between the 
broad Ly (3 width of 3130±155 kms -1 and the Ha value 
of 2800±300 kms -1 measured along the rim by Tuohy 
et al. (1982). However, our line width is over twice the 
1300±200 kms -1 width measured for the broad Ha line 
by Smith et al. (1991). The reason for the difference is 
that the spectroscopic slit used by Smith et al. (1991) 
was placed E-W across SNR 0519—69.0, intersecting the 
bright knot on the eastern rim of the remnant (clearly seen 
in Figure 2). The spectra of Smith et al. (1991) were 
only deep enough to reveal broad Ha emission from this 
bright knot, where the blast wave has encountered a rela- 
tively dense region of the ISM and slowed. On the other 
hand, Tuohy et al. (1982) observed the fainter, more uni- 
form portions of the SNR 0519—69.0 rim, where consider- 
ably less deceleration has occurred. As shown in Figure 2, 
the position and orientation of the LWRS aperture placed 
the brightened emission knots entirely outside the slit dur- 
ing the FUSE observation. Therefore, our FUV spectrum 
samples the less decelerated portions of SNR 0519—69.0 
and should give Ly (3 widths closer to the results of Tuohy 
et al. (1982). 

3.4. SNR 0548- 704 

There are no emission lines detected in the FUSE spec- 
trum of this object (Figure 6). The broad Ha width mea- 
sured by Smith et al. (1991) along the eastern rim 
of SNR 0548-70.4 is 760±140 kms -1 , corresponding to 
shock speeds of 700—950 kms -1 . Since we easily detected 
shocks of similar speed in our LWRS FUSE spectrum of 



DEM L 71 in both Ly (3 and O VI, it appears that the blast 
wave region covered by our LWRS slit (Figure 1) is propa- 
gating into a very low density, mostly ionized medium. In 
that case, pressure conservation would require very high 
shock speeds at this position, meaning that aside from be- 
ing faint, the lines are probably also very broad. 

Interestingly, our FUSE aperture includes several bright 
clumps of shocked interstellar material near the projected 
center of SNR0548-70.4 (Figure 1). The lack of C III 
and O VI emission from the knots in our FUSE spectrum 
suggests that they are excited by slow radiative shocks 
(^90 kms" 1 ; Hartigan, Raymond & Hartmann 1987). In 
their narrowband imagery of SNR 0548—70.4 Tuohy et al. 
(1982) detected all of these clumps in Ha, but only de- 
tected [O III] emission from the clumps lying outside the 
FUSE aperture. In addition, the Chandra ACIS S im- 
age of SNR 0548-70.4 (Hcndrick, Borkowski & Reynolds 
2003) docs not show any obvious X-ray emission from 
clumps within the FUSE aperture. These features are 
all characteristic of slow radiative shocks and suggest that 
dense, localized material has been overrun by the blast 
wave on the near and/or far side of SNR 0548—70.4. 

4. TEMPERATURE EQUILIBRATION AND KINEMATICS 

The FUSE spectra of the Balmer-dominated LMC 
SNRs are particularly valuable in that they feature promi- 
nent emission lines of hydrogen, carbon and oxygen, al- 
lowing us to directly measure the degree of ion-ion equili- 
bration in fast non-radiative shocks. The diagnostic infor- 
mation available to us includes (1) the relative widths of 
the C III A977, Ly (3 and O VI AA1032, 1038 lines, which 
reflect both the relative temperatures of these ions and 
the shock speed (assuming we know the contribution of 
bulk Doppler broadening to the line profiles), and (2) the 
Ly (3 I O VI flux ratio, which reflects the electron-proton 
equilibration, shock speed and preshock fraction of neutral 
hydrogen. 

In the case of SNR 0519—69.0, the good agreement be- 
tween our measured broad Ly (3 width and the Ha width 
measured by Tuohy et al. (1982) from the edge of the 
SNR suggests that the FUV profiles are primarily broad- 
ened by thermal motions, i.e., that the emission is domi- 
nated by the limb-brightened edges where the bulk motion 
along the line of sight is small. Since the LWRS slit also 
covers most of SNR 0509—67.5, we can expect the broad 
Ly (3 profile from this SNR to also be shaped by thermal 
broadening. Finally, we can expect our MDRS spectrum 
of the DEM L 71 blast wave to also exhibit lines domi- 
nated by thermal broadening, since the edge of the blast 
wave was targeted in this observation. In all these cases, 
we have used the FUV line widths as a direct measure of 
the ion temperatures and have utilized the non-radiative 
shock code of Ghavamian et al. (2001, 2002) to predict 
the emission line fluxes. 

In contrast to the spectra mentioned above, the inter- 
pretation of our LWRS line profiles for DEM L 71 is less 
straightforward. In this case most of the FUV emission 
is sampled from the projected center of the SNR, where 
emission from a range of viewing angles (and hence bulk 
velocities) is sampled in the FUSE spectrum. Further- 
more, the fraction of the DEM L 71 interior covered by the 
LWRS slit is large enough to make a simple plane-parallel 
approximation inadequate for describing the two sides of 
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the shell. In the next section we describe a more realistic 
model where we approximate the SNR as an expanding 
spherical shell. Having used these models to disentan- 
gle the bulk and thermal contributions to the line broad- 
ening, we then compute the proper shock speed(s) and 
T p /To ratio(s) required to match the FUV line profiles for 
DEM L 71. Finally, we model the flux ratios for DEM L 71 
using our numerical shock code and constrain the electron- 
ion equilibration and the preshock neutral fraction. 

4.1. NUMERICAL MODELS 

We begin by describing our plane parallel shock models, 
then describe the spherical shell model for DEM L 71. We 
use the former models to predict flux ratios for the broad 
Ly (3 and O VI emission lines, and use the latter models 
to match the shapes of these lines in the LWRS spectrum 
of DEM L 71. 

Before presenting our models, we note that Heng & 
McCray (2006; hereafter HM06) have recently presented 
new calculations of Balmcr-dominated shock structures. 
In their models HM06 found that the fast neutral distri- 
bution evolves in shape over multiple charge exchanges, 
with the overall broad component profile being the sum 
of multiple profiles (this effect was ignored in our earlier 
works). They find that for Vs ^2000kms _1 the fast neu- 
trals become ionized before they have undergone enough 
charge exchanges to fully sample the bulk-shifted velocity 
distribution of the incident hot protons. The resulting fast 
neutral distribution is shifted by a bulk velocity substan- 
tially less than | Vs , and is narrower than predicted from 
the models of Chevalier, Kirshner & Raymond (1980) and 
our own models (Ghavamian ct al. 2001, 2002). HM06 
find this effect to be far less important at the shock speeds 
relevant to DEM L 71, so the shell model results presented 
for that SNR in Section 5.1.1 should remain valid. How- 
ever, in the cases of SNR 0509-67.5 and SNR 0519-69.0 
where shock speeds are expected to exceed 2000 kms -1 , 
our interpretation of the FUV line profiles changes signifi- 
cantly if we adopt the HM06 mapping between shock speed 
and broad hydrogen line width. In the HM06 calculations 
the range of shock speeds implied by a given broad H I line 
width is significantly narrower in the limits of minimal to 
full electron-proton equilibration (for Vs J>2000 kms -1 ) 
than in our line profile calculations (Ghavamian et al. 
2001). We describe these results in more detail in Sections 
5.2 and 5.3. 

4.1.1. SHOCK MODEL 

The numerical code described by Ghavamian et al. 
(2001; 2002) computes the ionization and temperature 
structure behind a non-radiative shock with arbitrary 
electron-ion and ion-ion equilibration. The original pur- 
pose of the code was to compute the broad-to-narrow flux 
ratio of the Balmer Ha and H/3 lines, including the en- 
hancement of the narrow Balmer lines by Lyman line ra- 
diative transfer. We have now extended the ionization 
structure calculation of this code to include heavier ele- 
ments up to and including silicon. In addition to the elec- 
trons and protons, our code now includes Coulomb equi- 
libration of the He, C and O ion temperatures behind the 
shock. 

At the large shock velocities expected in SNR 0519—69.0 
and SNR0509-67.5 (V s ^4000 kms" 1 ), collisional ion- 



ization of neutral hydrogen by He + ions and alpha parti- 
cles becomes important (Laming et al. 1996), so we have 
included these processes in the code via rates computed 
with the atomic close-coupling cross sections of Toshima 
(1994). We have also updated the ionization balance to in- 
clude charge exchange between neutral hydrogen and both 
He + ions and alpha particles, again utilizing the cross sec- 
tions computed by Toshima (1994). All our rate calcu- 
lations assume Maxwellian particle distributions and are 
computed in the center of mass frame of the incident and 
target particles in the manner described by Ghavamian 
(1999). The other update included in this shock code is 
the calculation of emission line fluxes in broad Ly 7, C III 
A977, broad Ly (5 and O VI AA1032, 1038 out the front 
of the shock. We computed the collisional excitation of 
these lines by electron, proton, He + and alpha particle 
impact. The atomic data used in our calculation utilize 
fits to collision strengths for O VI excitation by electron 
impact from the compilation of Merts et al. (1980). To 
compute heavy ion excitation rates we have integrated the 
proton and alpha particle excitation cross sections of O VI 
from Laming et al. (1996), again under the assumption 
of Maxwellian particle distributions. We have also com- 
puted Ly (3 and Ly 7 rates from He + and alpha particle 
impact from the scaled cross sections of Janev (1996). 

Although all of the observed emission lines arise from 
resonant transitions, we have ignored the radiative transfer 
of the FUV photons within the shock in our calculations. 
The integrated resonance scattering optical depth perpen- 
dicular to the shock front in the broad Ly /3, Ly 7, C III 
and O VI lines is generally small, ^0.03 for a preshock 
density of 1 cm~ 3 , preshock hydrogen neutral fraction of 

0. 5 and shock speeds in the range 600—5000 kms -1 . In 
addition, the very large Dopplcr widths of the FUV lines 
(-600-3000 kms" 1 for the SNRs considered here) en- 
sure that most photons propagating perpendicular to the 
shock front escape before being absorbed. The prominent 
absorption features observed at the centers of the emis- 
sion lines are almost entirely due to LMC and Galactic 
halo absorption. 

Before describing the spherical shell model for DEM L 71 
we consider the question of whether the ion-ion equilibra- 
tion can be measured separately from the electron-proton 
equilibration in the LMC Balmer-dominated SNRs. At 
shock speeds greater than 3000 kms -1 , FUV line emission 
from collisional excitation by protons, He + ions and alpha 
particles becomes increasingly important. This implies 
that the flux ratios can be sensitive not just to the electron- 
proton equilibration, but the proton-helium ion equilibra- 
tion as well. Since there are no He lines in our FUSE 
spectra of the Balmer-dominated SNRs, we have no a pri- 
ori knowledge of the degree of proton-He ion equilibration 
in these objects. This raises the question of whether flux 
ratios uniquely determine the electron-proton equilibra- 
tion. The F(Ly /3)/Foy/(1032) ratio is most sensitive to 
Tue/Tp in shocks with minimal electron-proton equilibra- 
tion. For (T e /T p )o^0.05, the ratio F(Ly f3)/F y/(1032) 
for shocks in the range 3500—8000 kms -1 (appropri- 
ate for SNR 0509-67.5 and SNR 0519-69.0) increases by 
only 10% as Tn e /T p is varied from 4 (i.e., mH e /m p ) to 

1. Therefore, while the FUV line fluxes are sensitive to 
V s and (T e /T p )o in SNR 0509-67.5 and SNR 0519-69.0 
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(Vs <;3000 kms -1 ), they do not vary strongly with ion- 
ion equilibration over the range of shock speeds set by the 
width of the broad Ly (3 profile. 

4.1.2. SPHERICAL SHELL MODEL 

To predict the shapes of the FUV line profiles in our 
FUSE spectrum of DEM L 71 we have computed the 
emergent intensity in broad Ly (3 and O VI from two 
limb-brightened, hemispherical shells joined together in 
the plane of the sky. For each shell the shock speed, 
prcshock density and degree of temperature equilibration 
between the protons and oxygen ions are free parameters. 
From Figure 1 we estimate that approximately half the 
DEM L 71 shell is subtended by the LWRS slit. 

The first step in computing the line profile within the 
aperture is to compute the emission at a velocity v and an 
angle 9 from the center of a given shell via the well known 
radiative transfer relation 

I v (9) = S(l - e-^W) (1) 

where S is the source function of the shell and t is the 
optical depth of the shell at the azimuthal angle 9. Un- 
der the assumption of a Gaussian line profile, the optical 
depth for the front or back side of the shell is 

t v {9) = t exp(-(V T V B cos9f /b 2 ) sec 9 (2) 

where Vb = § Vs is the bulk velocity of the postshock gas, 
and b = (2/cT/M) 1 / 2 is the Doppler parameter of the 
emission line. The =F symbol is negative for the approach- 
ing shell and positive for the receding shell. The face-on 
optical depth to is evaluated at the center of the line pro- 
file, 

£*fj&L N (3) 

m e c b 

where N is the column density through the shock of the 
ion and are the oscillator strength and \j the rest 
wavelength of the ionic resonance line. In terms of these 
quantities, the source function S is 

In the limit To sec 9 <C 1, Equation (1) reduces to the usual 
expression for limb-brightened emission I ~ Iq sec 9, where 
Io is the face-on emission line intensity. If the shock speed, 
prcshock density and postshock equilibration are allowed 
to differ between the approaching and receding sides of 
the shell then S, Vb , b and To will also differ between 
the two sides. Although photons from the far side of the 
shell must traverse the near side of the shell as they exit 
the SNR, we assume that there is negligible absorption of 
these photons on the near side due to large velocity offsets. 
Despite the moderate blending of the Ly (3 and O VI lines 
in the LWRS spectrum of DEM L 71, the optical depth 
for absorption of highly Doppler shifted photons from one 
profile into another is negligible within the shock, so we 
neglect this effect in our models. 

The flux per velocity interval from each shell F(v) is 
obtained by inserting Equation (2) into Equation (1) for 
each of the two sides, then performing the integral 



F(v) oc / I v (9)cos9 d(cos9) (5) 
Jo 

for a sample of velocities over the line profile (we assume 
azimuthal symmetry). The total observed profile is then 
the summed fluxes from the front and back sides of the 
shell. 

In a spherical shell of radius R and thickness d, the 
maximum limb brightening is sec 9 max y // R/2d, where 
R is the radius of the SNR and d the thickness of the 
emitting layer. The larger the radius and the thinner 
the emitting layer, the greater the limb brightening. In 
our models we relate the upper limit on the angular inte- 
gral in Equation 5, 0, to the limb brightening factor via 
= a6 max , where 0< a < 1 is determined by the frac- 
tion of the SNR shell covered by the spectrograph slit. For 
a slit covering the full SNR, a = 1 . For a preshock density 
of 1 cm~ 3 , for example, do vi ~ (3-60) xl0~ 3 pc between 
shock speeds of 500 and 5000 kms" 1 . Taking R w 10 
pc for DEM L 71, the limb brightening factor lies in the 
range 9 <; sec 9 max <; 40, so that the integrated emission 
line profile F(v) is dominated by the contribution near the 
edge of the shell. 

To predict the line profile shapes we have used our nu- 
merical shock code to compute a grid of values for To over 
a wide range of shock speeds and equilibrations. When 
computing the line profile for a given set of shock param- 
eters, we interpolated within the pre-calculated grid to 
compute the appropriate value of t . For a given shock 
speed between 500 kms -1 and 5000 kms -1 the O VI 
column density behind the shock is relatively insensitive 
to the shock speed, but declines with increasing electron- 
proton equilibration, dropping by a factor of 3 between the 
limits of minimal and full electron-proton equilibration. 

For the range of shock parameters quoted above, the line 
center optical depths of broad Ly (3, O VI A1032 and C III 
A977 in our models can reach values ~0.5-l at the edge 
of the SNR shell. At these moderate optical depths the 
resonance scattering removes photons from the line cen- 
ter and leaves emission line profiles with slightly flattened, 
broadened peaks. In the case of the O VI doublet the 
F(A1032):F (A1038) ratio drops below the optically thin 
ratio of 2:1, an effect that has already been observed in 
non-radiative shocks in the Cygnus Loop (Long et al. 
1992; Sankrit & Blair 2002; Raymond et al. 2003). Al- 
though consideration of resonance scattering is important 
for interpreting the spectra of a localized filament in a 
nearby object such as the Cygnus Loop, the net effect on 
the spectra of the LMC Balmcr-dominated SNRs is rather 
small. In cases where the FUSE aperture includes edge- 
on emission from the SNR, enough O VI A1032 photons 
are contributed from positions just inside the edge of the 
shell to compensate for photons scattered out of the line 
of sight from positions at the edge of the shell. Our sim- 
ple shell models do not include these effects, so they may 
over-estimate the importance of resonance scattering on 
the FUV profiles of the LMC remnants. In comparing line 
fluxes from these objects with predictions from our plane 
parallel shock models we have neglected the effects of res- 
onance scattering. 

The line profiles calculated with the spherical shell 
model are flatter on top and exhibit slightly more ex- 
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tended wings than Gaussian profiles. However, given the 
strong interstellar absorption near the centers of the FUV 
lines and the limited signal-to-noise of most FUV line 
profiles, both the single Gaussian profile and the spher- 
ical shell model adequately match the line shapes in the 
SNR 0519—69.0 and SNR 0509—67.5 spectra. 

5. MODEL RESULTS 

In previous modeling of Balmcr-dominated SNRs 
(Ghavamian ct al. 2001; 2002), we estimated the degree 
of electron-proton temperature equilibration by first mea- 
suring the width of the broad Ha component from edge- 
on shocks to constrain the range of shock speeds, then 
modeling the broad to narrow flux ratio to constrain the 
combination of shock speed and (T e /T p )o consistent with 
all the data. Unfortunately, the narrow Ly (3 emission is 
completely absorbed by the ISM in our FUSE spectra, as 
is the core of the broad Ly (3 line. In addition, the widths 
of the FUV emission lines in all our FUSE spectra reflect 
both thermal and bulk Doppler broadening. Without tak- 
ing into account the latter contribution to the line profile, 
the shock velocities tend to be overpredicted. 

In the ensuing analysis we take the following approach 
to constraining the shock kinematics: (1) We use existing 
measurements of the broad component Ha width (where 
available) from the rim of each SNR to constrain the range 
of shock speeds. The line widths from these measure- 
ments should primarily reflect thermal broadening and 
give the best constraint on the appropriate range of Vs ■ 
(2) Using the range of constrained shock speeds from the 
plane parallel models, we calculate fluxes in broad Ly (3 
and O VI. Comparing these model fluxes with the ob- 
served values, we attempt to further refine our estimates of 
Vs and (T e /T p )o, while also obtaining constraints on the 
preshock neutral fraction of hydrogen. (3) Again utilizing 
the shock speed constraints from the plane parallel mod- 
els, we use our spherical shell models to predict the shapes 
of the integrated Ly (3 (and when possible, O VI) profiles. 
Comparing these predictions with observed profiles in the 
LWRS spectra, we gauge the contribution of bulk Doppler 
broadening to the profile width. By subtracting the bulk 
Doppler component in quadrature from the total profile 
width in Table 2, we obtain limits on T p and To, thereby 
constraining {To/T p )q. 

In the following we compare the observed F(Ly 
P)/Fo v/(1032) with predictions from our plane paral- 
lel shock models to constrain (T e /T p ) for the DEM L 71 
blast wave, SNR 0509-67.5 and SNR 0519-69.0. In all 
calculations we adopt the LMC abundances of Russell & 
Dopita (1992) (in this case by 12 + log(0/H) = 8.35). 

5.1. DEM L 71 

5.1.1. SOUTH RIM OF DEML 71 

The FWHM of the broad Ly (3 line in our MDRS spec- 
trum is 1365±75 kms" 1 , nearly 40% larger than the Ha 
width of 805i"£ km s 1 measured by Ghavamian et al. 
(2003) from Rutgers Fabry-Perot (RFP) observations of 
DEML 71. There are several potential explanations for 
this discrepancy First, there was a lack of extensive spec- 
tral coverage on either side of the Ha profiles in the RFP 
spectra of Ghavamian et al. (2003). Since the Ha pro- 
files were not sampled beyond ±800 kms -1 , fitting a line 



with wings extending well beyond that range may have 
resulted in an anomalously low line width. On the other 
hand it is possible that contrary to our expectation, the 
MDRS observation of DEM L 71 did not completely isolate 
the limb emission. In recent optical spectra of DEM L 71 
acquired with the Magellan 6.5m telescope (C. Rakowski 
2006; private communication) several broad Ha compo- 
nents are observed along portions of the eastern rim. The 
two components are of comparable width (^600 kms -1 ), 
but their line centroids are shifted in opposite directions 
in velocity. This indicates that even near the limb of the 
SNR the slightly off-edge emission can be bright enough to 
significantly affect the observed broad component profile. 
The multiple nested filaments observed at the southern 
limb of DEML 71 (Figure 2) suggest multiple tangencies 
between the blast wave and line of sight, indicating that 
the emission between these tangencies is likely to include 
partially face-on emission. Unlike the RFP spectral ex- 
traction, which cleanly isolated the Ha emission from the 
southern rim of DEML 71, the MDRS slit is wide enough 
to include some of this off-tangent emission, so the MDRS 
spectrum is more likely to be affected by this problem. If 
this explanation is correct, then the line widths measured 
in our MDRS spectra would overpredict the shock velocity. 

With the above caveats in mind, we modeled the 
DEML 71 MDRS spectrum by first using the Ha line 
width of 8051^ kms" 1 from Ghavamian et al. (2003) 
to set the range of shock speeds. Since the shock speed 
is likely to remain relatively constant near the rim of 
DEML 71, we assume that the FUV line ratios reflect ex- 
citation from a single shock, even though the lines exhibit 
some bulk broadening. Under these assumptions, we then 
compared the predicted F(Ly (3)/Fo v/(1032) ratios with 
the observed values. 

Using the relationship between hydrogen line FWHM 
and shock speed calculated by Ghavamian et al. (2001), 
the implied range of shock speeds from the Ha broad 
component width (Ghavamian et al. 2003) is 775—1005 
kms -1 over the range (T e /T p ) — m e /m p to 1. We calcu- 
lated a grid of F(Ly /3)/Foy/(1032) values over this range 
of shock speeds, combining each (T e /T p )o with the corre- 
sponding velocity needed to match the width of the broad 
Ly (3 line. A plot of our results is shown in the top panel 
of Figure 8. 

For the range of shock speeds here, the electrons in the 
H I neutral zone behind the shock are close to the opti- 
mum temperature (^10 6 K) for exciting Ly (3 close to the 
shock front, producing a peak in the F(Ly /3)/Foy/(1032) 
curve near (T e /T p ) =0. At higher equilibrations and 
shock speeds, however, the charge exchange rate begins 
to dominate over the collisional ionization rate, resulting 
in higher overall density of fast neutral hydrogen behind 
the shock. The broad Ly (3 flux increases accordingly, pro- 
ducing the rise in F(Ly (3)/Fo y/(1032) seen in Figure 8 
at high equilibrations. 

Unfortunately, due to the linear dependence of F(Ly 
(3)/Fo y/(1032) on preshock neutral fraction, the double 
valued behavior of the predicted F(Ly /?)/Foy/(1032) at 
low and high (T e /T p )o, and the uncertainty in the mea- 
sured flux ratio, we are unable to obtain a unique com- 
bination of (T e /T p )o and Vs from the plot. Although 
we cannot unambiguously determine these two quantities, 
we can at least constrain the range of preshock hydrogen 
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neutral fractions in the southern rim of DEM L 71. If we 
require at least one point of intersection between the range 
of measured and predicted F(Ly (3) /Fo y/(1032), the max- 
imum and minimum allowed preshock neutral fractions are 
40% and 20%, respectively. 

Assuming a uniform shock speed along the southern 
rim of DEM L 71, we have attempted to estimate the de- 
gree of proton-oxygen equilibration at the shock transition. 
At the range of shock speeds considered here (800—1000 
kins" 1 ), charge exchange behind the shock produces a fast 
neutral distribution with a velocity width nearly equal to 
the width of the proton distribution. Therefore, the pro- 
ton temperature is just the temperature obtained from 
the broad component Ha width (805 kms -1 ). If we 
assume that the difference between this value and our 
measured Ly (3 width is due to the effect of bulk mo- 
tion of emitting material within the FUSE aperture, then 
we can estimate the Doppler broadening as ~ (1365 2 
- 805 2 ) 1 / 2 = 1100 kins" 1 . This value, however, is 
larger than the measured O VI line width of 935 kms -1 . 
Therefore, To is not well constrained from our fit to the 
O VI profile, and we can only estimate an upper limit on 
{To /T p )ounder the assumption that the O VI width repre- 
sents pure thermal broadening. Under these assumptions, 

(To/T p )o = 4f*> giving (T /T p )o £ 8. Note that the 

oxygen temperature derived from the O VI line width rep- 
resents an average value weighted by the postshock O VI 
emissivity. The O VI zone behind the shock is extended 
enough for Coulomb collisions to have transferred some en- 
ergy from the oxygen ions to the lighter ions. Comparing 
our measured values of (To/T p ) with predictions from 
shock models, we find the emission-weighted ratio should 
be ~10% smaller than the immediate postshock value, so 
that (To/T p )o <J 9 is a better estimate of the temperature 
ratio. Clearly this is a rough limit, but it suggests that 
a moderate amount of proton-oxygen temperature equili- 
bration has taken place at the shock transition. 

5.1.2. LWRS SPECTRUM OF DEM L 71 

As described earlier, the LWRS slit in the DEM L 71 ob- 
servation is filled mostly with emission from the interior 
of DEM L 71. Therefore, we expect each of the emission 
lines in Figure 3 to consist of two components, one red- 
shifted and one blue-shifted from the rest frame of the 
SNR. Utilizing our spherical shell model for DEM L 71 we 
have attempted to match the shapes of the broad Ly (3 and 
O VI A 1032 lines. The faintness of the C III A977 feature 
prevents us from obtaining meaningful constraints on its 
shape, so we have not attempted to model the C III A977 
profile. 

We focused first on separating the thermal and geomet- 
ric contributions to the width of the broad Ly (3 line. We 
attempted to match the profile shape with our spherical 
shell models by varying Vs and (T e /T p ) together for 
each of the two sides of the shell contained in the aper- 
ture. The amplitudes of the frontside and backside pro- 
files arc not allowed to vary independently, but are tied 
together via the requirement of ram pressure conservation: 
nbvl — n f v 2 j . We do not model the LMC and Galactic 
Halo absorption in these models, but focus upon matching 
the wings of the observed Ly (3 profile. 

Our best estimate of the Ly (3 profile shape from the 



shell models is shown in Figure 9. As a first estimate, we 
assumed the LWRS slit was aimed at the center of the SNR 
and covered half the projected area of that shell (approxi- 
mately correct for our observational setup). Interestingly, 
the line profiles predicted by this model (dashed line in 
Figure 9) overpredict the emission near the center of the 
observed profile, and give only slightly better agreement 
in the far wings of the profile. Varying the amplitude or 
width of the model profile for any individual component 
(approaching or receding side of the shell) can produce bet- 
ter agreement toward the center of the broad Ly (3 line, 
but the agreement then worsens in the wings of the pro- 
file. Since the profiles from the front and back sides of the 
shell are centered on the bulk velocities (=§Vg ) of these 
two components, the blueshifted and rcdshifted peaks are 
spread far apart, with too little emission predicted at the 
center of the profile. We found that the only solution to 
this problem was to include emission from edge-on shocks 
in the profile calculation (i.e., by setting a = 1 in Equa- 
tion 5), despite the fact that the LWRS observation does 
not appear to include filaments near the edge of DEM L 71 
(Figure 1). Under that assumption the wings of the cal- 
culated profile draw inward, with emission from zero ra- 
dial velocity filling in the center of the broad Ly (3 line. 
The necessary inclusion of an edge-on component in our 
models indicates that some of the O VI emission may be 
distributed differently from the Ha emission (particularly 
if parts of the DEM L 71 blast wave are propagating into 
fully ionized gas). 

While the addition of the edge-on component signifi- 
cantly improves the ability of our models to match the 
observed profile, asymmetric flux contributions from the 
front and back sides of the shell are also required to match 
the residual redshift of the broad Ly (3 centroid. To repro- 
duce the +235 kms -1 redshift of the broad Ly (3 feature 
from the LMC rest velocity (Table 2), the central peak 
of the Ly (3 profile from the far side of the shell must be 
higher that of the near side. Keeping this constraint in 
mind, we adjusted the shock velocity, total preshock den- 
sity and electron-proton equilibration for the two sides of 
the shell while conserving ram pressure as described ear- 
lier. The final result of our Ly (3 profile calculation for 
DEM L 71 is marked with the solid line in Figure 9. 

Using the calculated relationship between the velocity 
width of the proton distribution and the width of the H I 
profile transmitted by charge exchange (Chevalier, Kirsh- 
ner & Raymond 1980; Ghavamian et al. 2001), the best 
combination of Ly (3 profile widths from our spherical shell 
model are 920 km s -1 for the approaching side of the shell 
and 585 kms^ 1 FWHM for the receding side of the shell. 
The corresponding combination of shock parameters (Vs , 
n H o, (T e /T p ) ) are (900, 0.5, 0.05) for the approaching 
side of the shell and (650, 1.0, 0.6) for the receding side of 
the shell. Note that the absolute quoted values of n H o are 
not the meaningful quantities here, but rather their ratio, 
as required by ram pressure conservation. 

Although there is some degeneracy between the vari- 
ous shock parameters estimated above, they are fully con- 
sistent with prior measurements of the broad component 
Ha width, shock speed and equilibration estimated for 
DEM L 71 by Smith ct al. (1991), Ghavamian ct al. 
(2003) and Rakowski et al. (2003). The analyses in 
these papers indicated that the blast wave of DEM L 71 
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has decelerated to ^600 km s in the northeastern corner 
of the SNR, precisely where the FUSE LWRS slit grazes 
the edge. On the other hand, the RFP Ha spectra and 
Chandra X-ray analyses of Ghavamian et al. (2003) and 
Rakowski ct al. (2003) indicated a nearly 50% variation 
in shock speed around the rim of of DEM L 71, with the 
highest shock speeds (~800-1000 kms -1 ) measured along 
the southern rim. The results of our broad Ly (3 profile 
models indicate that the far side of DEM L 71 may be 
propagating into a higher density medium than the near 
side, with a correspondingly lower shock speed on the far 
side of the SNR. 

After applying our shell models to the broad Ly (3 pro- 
file, we modeled the O VI A1032 profile. As a first guess 
we set the shock speeds of the front and back sides of 
the shell equal to the values derived from the broad Ly (3 
profile model, with the O VI line width a free parameter. 
Our best estimates for the O VI line widths are 475 km s _1 
FWHM for the receding shell and 900 kms" 1 FWHM for 
the approaching shell. The centroid of the observed O VI 
A1032 feature is shifted by the same velocity as the cen- 
troid of the Ly (3 feature (Table 2), and can be matched 
with a 2:1 ratio of preshock density between the back and 
front sides of the shell. Taken together, the Ly f3 and O VI 
A1032 results suggest that while the front and back sides 
of DEM L 71 are propagating into an ISM with a 2:1 ratio 
in total density, the hydrogen neutral fractions are similar 
on both sides of the shell. 

From our above analysis, we can estimate (To/T p )q 
for the approaching and receding sides of the shell in 
DEM L 71. Given V s =650 kms" 1 , (T e /T p ) = 0.6 and 
an estimated O VI A1032 width of 475 kms -1 on the 
back side, we obtain (To/T p )q « 6.5. On the near side, 
for shock speed of 900 kms" 1 , (T e /T p ) = 0.05, we ob- 
tain (To /T p )o~ 16. Although these estimates are sub- 
ject to the uncertainty of deconvolving the thermal broad- 
ening from the bulk Doppler broadening in our spectra, 
both the magnitude and declining trend of (To/T p )o with 
shock speed are consistent with results established for 
the Balmer-dominated shocks in other SNRs such as the 
Cygnus Loop (Raymond et al. 2003) and SN 1006 (Lam- 
ing et al. 1996; Korreck et al. 2004). 

5.2. SNR 0509-67.5 

The sole emission line positively detected in the LWRS 
spectrum of this Balmer-dominated SNR is broad Ly f3 
(Figure 4). The fitted FWHM (3710±400 kms" 1 ) is 
the largest hydrogen line width measured in a Balmer- 
dominated SNR to date. The centroid of the broad Ly 
(3 line is substantially blueshifted (—540 kms" 1 ) relative 
to LMC frame, suggesting that the FUV emission from 
SNR 0509—67.5 primarily arises from the near side of the 
SNR shell. It is likely that the broad Ly [3 profile is affected 
by bulk Doppler motions. However, unlike the other three 
Balmer-dominated SNRs, there has been no prior detec- 
tion of broad ionic line emission from SNR 0509—67.5, so 
we have no clear constraint on the appropriate range of 
shock speeds for modeling the profile in Figure 4. 

Since we do not have enough information for disentan- 
gling the bulk Doppler contribution of broad Ly (3 from the 
thermal Doppler contribution, we proceed with a plane 
parallel shock analysis under the crude assumption that 
the broad Ly (3 width in SNR 0509-67.5 reflects purely 



thermal broadening. Using the calculations of HM06, 
the measured FWHM corresponds to shock speeds of 
5200-6300 km s" 1 over the range m e /m p < (T e /T p ) < 1 
(for comparison, the FWHM calculations of Ghavamian 
et al. (2001, 2002) give a range of 5610-7780 kms" 1 ). 

Although the O VI lines are not detected in the FUSE 
spectrum of SNR 0509—67.5, we can still place a lower 
limit on both F(Ly /3)/F O y/(1032) and the widths of the 
O VI lines. First, we estimate an upper limit on the O VI 
A1032 flux. The upper limit on the flux per wavelength 
interval from the FUSE spectrum of SNR 0509—67.5 is ~ 
2xl0~ 15 ergs cm~ 2 s _1 A -1 at the center of the O VI 
A1032 line. To convert this number to a flux, we require 
a limit on the O VI line width. For each shock speed in 
the range quoted earlier, the FWHM of the O VI line is 
related to the oxygen temperature via Vfwhm{0 VI) = 
V8 2 \JkTo I mo , so that full proton-oxygen equilibra- 
tion would produce an O VI line width j as broad as the 
proton velocity distribution, while no equilibration would 
yield equal widths. Over the full range of allowed equili- 
brations, the predicted O VI line width ranges from 1245 
kms" 1 (for(T e /Tp) = (To/T p )o = l)to5300kms- 1 (for 
(T e /T p )o = 1 and (To/T p )q = 16). kms -1 (equilibrated). 
The corresponding limit on the O VI A 1032 flux is then 
F(1032) = (0.9-3.9) xl0~ 14 ergs enr 2 s" 1 in the limits of 
full and minimal proton-oxygen equilibration respectively. 

In Figure 8 we present the predicted F(Ly 
/3)/F O y/(1032) ratio for SNR 0509-67.5 for the range 
of shock speeds implied by the width of the broad Ly (3 
profile. The lower limit on the F(Ly /3)/F oy/ (1032) ratio 
corresponds to the largest O VI line width (^5300 kms -1 ) 
allowed from our analysis above. For a given broad Ly (3 
width and corresponding combination of shock parameters 
((T e /T p ) , Vs ), the O VI emission rises more rapidly with 
increasing (T e /T p )o (and shock speed) than the broad Ly 
(3 emission, causing the modest decline in the predicted 
F(Ly /3)/F O y/(1032) (Figure 8). Our models are shown 
for two limiting preshock neutral fractions: fn<> = 0.4 and 
0.8. The lower neutral fraction is the smallest value allow- 
ing an overlap between our model curves and the observed 
range. The higher fraction is the largest value likely for 
the warm ISM, given the minimum photoionization ex- 
pected from both the SNR and the ambient UV and soft 
X-ray backgrounds. 

We can draw several conclusions from the middle panel 
of Figure 8. First, efficient proton-oxygen ion equilibration 
is ruled out. In such a case, the narrow O VI line width (a 
minimum of 1245 km s" 1 FWHM for (T /T p ) = 1) would 
require a very low total O VI A1032 flux to stay consistent 
with the non-detection of O VI in the LWRS spectrum. 
This would raise the observed F(Ly /3)/F O y/(1032) ratio 
(i.e., the lower edge of the shaded region) to <~ 9, which 
cannot be matched by any model regardless of equilibra- 
tion or preshock neutral fraction. The prediction from 
our models is that F(Ly [3)/¥o y/(1032) ^4.5, a constraint 
which is satisfied for V FW hm(0 VI) ^ 2400 kms -1 . This 
places the proton-oxygen ion equilibration in the range 
4^ To/Tp < 16 (the assumed upper limit corresponds to 
the oxygen to proton mass ratio, though (To/T p ) > 16 
has been observed on occasion in solar wind shocks). We 
cannot place a strong constraint on (T e /T p )o from Fig- 
ure 8, but we can deduce that the preshock gas should be 
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significantly neutral, 0.4 <J /#o <^0.8. 

We can use our derived shock velocity for 
SNR 0509—67.5 to estimate its age. Assuming pure ther- 
mal broadening for the detected Ly line and taking into 
account the measurement uncertainty in the broad Ly 
width, the predicted shock speed is 3600—7100 kms -1 
over the range m e /m p < (T e /T p ) < 1. We compute the 
age r by assuming a radius of 3.6 pc (Warren & Hughes 
2004) and using t = 5 y, where S is the expansion pa- 
rameter. Between the limits of free expansion and the 
asymptotic Sedov- Taylor limit, S ranges from around 1.0 
down to 0.4. Our derived shock speeds are significantly 
lower than the 10,000-20,000 kms -1 expansion veloc- 
ities predicted for Type la SNe, indicating that some 
deceleration of the blast wave has occurred. Assuming 
standard Type la parameters and a constant ISM density, 
Dwarkadas & Chevalier (1998) predicted that S declines 
to around 0.6 during the first few hundred years after the 
SN explosion. Using this expansion parameter, we find 
t w 295-585 years for SNR 0509-67.5. This value com- 
pares favorably with the light echo age of 400±120 years 
found for SNR by Rest et al. 2005). Given the neglect 
of bulk Dopplcr broadening in our analysis, we are at the 
very least encouraged by the broad agreement between 
our derived age for SNR 0509—67.5 and the independent 
age estimate of Rest ct al. (2005). 

5.3. SNR 0519-69.0 
5.3.1. PLANE PARALLEL MODEL RESULTS 

The broad component Ha width measured along the 
rim of SNR 0519-69.0 by Tuohy et al. (1982) was 
2800±300 kms -1 , comparable to the value measured in 
the Galactic Balmer-dominated SNR SN 1006 (Ghavamian 
et al. 2002; Korreck et al. 2004). From the line 
profile calculations of HM06 this width corresponds to 
shock speeds in the range 3200—3900 kms -1 in the lim- 
its m e /m p < (T e /T p )o < 1 (again these are considerably 
narrower values than the range of 4115—5930 km s _1 pre- 
dicted by our own models). 

The predicted F(Ly /3)/F O y/(1032) ratios for 
SNR 0519— 69.0are shown in the bottom panel of Fig- 
ure 8. It is clear that in this case the flux ratios do not 
vary strongly with electron-proton equilibration, partly a 
result of the narrow range of shock speeds implied by the 
width of the broad Ly line. The F(Ly /3)/F OVJ (1032) 
curves are shown for two limiting preshock neutral frac- 
tions: the lower fraction (fgo = 0.4) being the smallest 
value allowing overlap between our model curves and the 
observed range, with the higher fraction (f H o = 0.8) being 
the largest value likely for the warm ISM surrounding the 
SNR. We cannot place a strong constraint on (T e /T p )o 
from Figure 8. However, the preshock neutral fraction is 
well constrained, 0.4 ^ f H o ^0.5. 

Despite the considerable uncertainty (~40%; Table 2) 
in the O VI line width, we can still place reasonable limits 
on the degree proton-oxygen ion equilibration in the blast 
wave of SNR 0519-69.0. The FWHM of the broad Ly 
line is 3130±155 kms" 1 , while that of the O VI A1032 
is 4975±1830 kms" 1 . The contribution to the FUV pro- 
files from bulk Doppler broadening, Vb, can be estimated 
by subtracting the Tuohy et al. (1982) Ha width from 
the observed Ly width, V B « V3130 2 - 2800 2 = 1400 



km s 1 . Assuming the same bulk broadening for the O VI 
line, the thermal width of the oxygen ion distribution is 
then 4800±1750 kms -1 . Since the O VI line width cannot 
significantly exceed that of broad Ly 0, our result can only 
be consistent with equal line widths, or (To/T p )q ps 16. 

The age obtained by Rest et al. (2005) from their light 
echo study of SNR 0519-69.0 is 600±200 years. Taking 
into account the measurement uncertainty in the broad Ha 
width, the full range of shock speeds for SNR 0519—69.0 
between the limits of minimal and full electron-proton 
equilibration is 2600—4500 kms -1 . Taking S w0.6, this 
results in an age of 520—900 yrs, in good overall agreement 
with the age predicted by the light echo study. 

5.3.2. SPHERICAL SHELL MODEL RESULTS 

Using the above constraints on the shock speed and equi- 
libration, we have computed the shape of the broad Ly 
line produced by a spherical shell. Using Figure 7 of 
HM06 we find that for an Ha line width of 2800 kms" 1 
(the edge-on shock value of Tuohy et al. 1982) the 
predicted bulk velocity of the fast neutrals is ~ jVs, 
rather than the usually assumed value of | Vs . In that 
case 2000^ V B ^OOOkms- 1 for SNR 0519-69.0 in the 
range m e /m p < (T e /T p )o < 1. In addition, since the cen- 
troids of the FUV emission lines in SNR 0519-69.0 (Ta- 
ble 2) are consistent with the LMC rest velocity, we assume 
equal shock speeds and preshock densities for the front 
and back hemispheres. Using all of these parameters, we 
were able to match the observed broad Ly profile with 
a thermal broadening FWHM of 2500 kms -1 , consistent 
with the lower bound on the Ha width of Tuohy et al. 
(1982). The corresponding shock speed is approximately 
2600 kms -1 from Figure 7 of HM06, with a correspond- 
ingly low electron-proton equilibration ((T e /T p )o <^0.2). 
This result is consistent with shock speeds at the lower end 
of the range predicted in the previous section. It worsens 
the disagreement between the light echo age of Rest et al. 
(2005) and our dynamical age estimate, and may provide 
further evidence of blast wave deceleration. 

6. THE EFFECT OF ASSUMED ABUNDANCES 

It is clear that our interpretations above are significantly 
affected by our choice of oxygen abundance via the lin- 
ear dependence of the F(Ly /3)/Foy/(1032) ratio on the 
O/H abundance ratio. Hughes, Hayashi & Koyama (1998) 
obtained a lower average oxygen abundance than Russell 
& Dopita (1992) from their analyses of ASCA spectra of 
LMC SNRs: 8.21±0.07 dex. It is instructive to consider 
the effect of such lower O/H ratios on the comparison of 
our shock models with the observations. If ISM oxygen 
abundances around all four Balmer-dominated SNRs were 
set to the values derived by Hughes, Hayashi & Koyama 
(1998), the modeled O VI A1032 flux would decrease by 
28%, raising the predicted F(Ly 0)/F o y/(1032)ratios in 
Figure 8 and lowering the range of acceptable preshock 
neutral fractions. In the case of DEM L 71, the pre- 
dicted range of acceptable preshock neutral fractions be- 
comes 15%-30%, while a fraction above 30% is allowed for 
SNR 0509-67.5. The neutral fraction for SNR 0519-69.0 
is most tightly constrained, with a value close to 30% re- 
quired to match the new F(Ly 0)/Fq y/(1032)ratio. 
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7. REDDENING CONSTRAINTS FROM THE FUV SPECTRA 
OF DEM L 71 AND SNR 0519-69.0 

Aside from the kinematic diagnostics available in the 
FUV spectra of the Balmer-dominated SNRs, the ratio of 
Ly (3 to Ly 7 flux is sensitive to the interstellar redden- 
ing, providing further useful information on these objects. 
Over the range of Vs and (T e /T p ) considered in this pa- 
per, the impacting electrons and ions have kinetic energies 
greatly exceeding the energy difference between the n = 3 
and n = 4 levels of hydrogen, and the ratio of these two 
lines reduces to the ratio of their collision strengths. In 
addition, the Ly (3 and Ly 7 lines are optically thin due 
to their large Dopplcr widths and the low density of the 
postshock gas. This allows us to ignore radiative transfer 
of Ly (3 and Ly 7 photons when computing their flux ratio. 
The ratio is then nearly constant, ranging from approxi- 
mately 3.3 for minimal equilibration to 3.1 for maximal 
equilibration. 

In the LWRS spectra of DEM L 71 and SNR 0519-69.0, 
where emission from both lines has been detected, the 
Ly (3 I Ly 7 ratios are 3.7±0.4 and 3.4±0.5, respectively. 
Comparing these numbers with the predicted values, these 
numbers are consistent with low to minimal reddening. 
Koornneef & Code (1981) noted that the slopes of the 
Galactic and LMC extinction curves are very similar at 
A < 1400 A. This implies that for the ratio of two line 
fluxes in this wavelength range, the correction will be 
nearly identical whether we use an LMC or Galactic ex- 
tinction law. Indeed, the range of extinctions we derive for 
DEM L 71 and SNR 0519—69.0 using a Galactic redden- 
ing law (Cardelli, Clayton & Mathis 1989) are very simi- 
lar to the range derived for the non-30 Doradus regions 
of the LMC (Fitzpatrick 1985). The color excess pre- 
dicted for DEM L 71 is then E(B - V) = 0.05-0.13, while 
for SNR 0519-69.0 we find E(B-V)<0.11. 

8. DISCUSSION AND CONCLUDING REMARKS 

We have presented the first far ultraviolet spectra of the 
four known Balmer-dominated SNRs in the LMC. These 
objects - DEM L 71, SNR 0509-67.5, SNR 0519-69.0 and 
SNR 0548—70.4, are the non-radiative (adiabatic) rem- 
nants of Type la SN explosions. They feature optical spec- 
tra dominated by Balmer line collisional excitation of am- 
bient H I gas overrun by the blast wave. Our FUSE obser- 
vations resulted in the detection of three of the four SNRs, 
DEM L 71, SNR 0509-67.5 and SNR 0519-69.0. The de- 
tected emission lines - C III A977, O VI AA1032, 1038, 
Ly (3 and Ly 7 — are very broad (600 — 3700 kms -1 ) 
and are generated in the thin ionization zones behind the 
non-radiative shocks. In the case of SNR 0509— 67.5, the 
only detected feature is faint, broad Ly (3 (FWHM of 3710 
kms" 1 ); however, this is the first and only broad emis- 
sion line detected in this SNR since its discovery by Tuohy 
et al. (1982). 

The relative fluxes of the FUV emission lines are sen- 
sitive to the shock velocity, preshock neutral fraction and 
the degree of electron-proton equilibration at the shock 
front. On the other hand, the widths of the FUV emis- 
sion lines in Balmer-dominated SNRs are proportional to 
the relative temperatures of the heavy ions, making them 
useful diagnostic probes of the degree of ion-ion temper- 
ature equilibration in collisionlcss shocks. These proper- 



ties make FUV spectra of these SNRs valuable tools for 
probing the collisionless heating in high Mach number col- 
lisionless shocks. 

Comparing the F(Ly (3)/Y y/(1032) flux ratios with 
predictions from plane parallel and spherical shell mod- 
els of non-radiative shocks, we find that the observed 
flux ratios are consistent to within the errors with low 
to minimal electron-proton equilibration in shocks faster 
than 2500 kms" 1 (observed in SNR 0509-67.5 and 
SNR 0519—69.0). On the other hand, spectra from the 
interior and limb of DEM L 71 (V s -600-1000 kms" 1 ) 
show evidence for small to moderate electron-proton tem- 
perature equilibration. Our results are generally consistent 
with the inverse trend between electron-proton equilibra- 
tion and Mach number derived from optical (Ghavamian 
et al. 2001, 2002; Rakowski 2005, Ghavamian, Laming & 
Rakowski 2007) and joint optical and X-ray (Ghavamian, 
Laming & Rakowski 2007; Rakowski, Ghavamian & 
Hughes 2003) observations of Balmer-dominated shocks in 
SNRs. 

We have estimated (To/T p ) at the shock transition by 
comparing the broad Ly (3 and O VI emission line widths 
in our FUSE spectra of the Balmer-dominated LMC rem- 
nants. Over the range of shock speeds sampled by our 
FUSE observations (~600^V S ^6000 kms" 1 ) our ob- 
servations indicate that (To/T p ) increases from <~ 7 for 
shock speeds near 600 kms -1 to ~ 16 (mass-proportional 
heating) for shock speeds exceeding 1000 kms -1 . A sum- 
mary of our modeling results is shown in Table 3. From 
the shock velocities listed in Table 3 we obtain ages of 
295-585 years for SNR 0509-67.5 and 520-900 years for 
SNR 0519—69.0, respectively, in good agreement with the 
ages estimated for these SNRs by SN light echo studies 
(Rest et al. 2005). 

The main source of uncertainty in our estimates of 
(To/T p )q is our uncertainty in the O VI emission line 
widths. In the case of SNR 0519—69.0, the faintness of the 
O VI lines and their heavy blending with broad Ly (3 makes 
a precise measurement of their widths difficult. In the case 
of SNR 0509-67.5, the O VI emission falls below the de- 
tection threshold, allowing only an indirect estimate of the 
O VI line widths. Although line blending is far less severe 
in our FUV spectrum of the southern rim of DEM L 71 all 
the emission lines are intrinsically very faint, again mak- 
ing the measurement of the O VI line widths (and hence 
(To/T p )q) uncertain. Line blending is also less severe in 
our FUSE spectrum of the interior of DEM L 71. However, 
disentangling the frontside from backside emission in this 
case is the dominant source of uncertainty in measuring 
(T o /T p ) . 

Modulo the bulk Doppler broadening, our results are 
consistent with an inverse relationship between the de- 
gree of ion-ion equilibration and shock speed, with slower 
shocks exhibiting more prompt equilibration at the shock 
front than faster shocks. This trend is similar to that 
found for the electron-proton temperature equilibration, 
suggesting that the collisionlcss processes responsible for 
prompt electron-proton equilibration at low shock speeds 
may also cause prompt ion-ion equilibration at low shock 
speeds. A stronger verification of this trend will require 
further measurements of (To/T p )q in other non-radiative 
shocks at intermediate velocities (500—2500 kms -1 ). Ob- 
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taining these observations will remain the subject of future 
work. 
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Table 1 
List of FUSE Observations 



Object 


Program ID 


aj20oo 


(5.72000 


Aperture 


Exp. (ks) a 


SNR 0505-67.9 (DEM L 71) 


P21403 


05 h 05 m 42!7 


-67° 52' 38'.'0 


LWRS 


27.7 (12.7) 


SNR 0505-67.9 (DEM L 71 SW) 


C072(03-04),C07206 


05 h 05 m 39!7 


-67° 53' 14'/4 


MDRS 


100.2 (44.6) 


SNR 0509-67.5 


P21402 


05 h 09 m 31!9 


-67° 31' 17'.'2 


LWRS 


29.6 (9.7) 


SNR 0519-69.0 


P21401 


05 h 19 m 33!8 


-69° 02' 09'.'7 


LWRS 


27.4 (12.7) 


SNR 0548-70.4 


P21404 


O5 h 47 m 50?l 


-70° 24' 52'/ 1 


LWRS 


29.5 (10.4) 



a Total integration time for each object. Orbital night time fractions arc listed in parentheses. 



Table 2 

Measured Parameters for LMC Balmer-Dominated SNR Spectra 1 



Object 



Broad Ly (3 O VI A1032 O VI A1038 b Broad Ly 7 C C III A977 



SNR 0505 67.9 (DEM L 71) 

P214 (LWRS) 

•Vent [A] 

Vccnt [km S" 1 ] 
Vfwhm [km s _1 ] 
Flux [10~ 14 ergs cm" 
C072 (MDRS) 

Acent [A] 

V C cnt [km s" 1 ] 
Vfwhm [km s _1 ] 
Flux [10 -14 ergs cm" 





1027.1±0.05 


1033.2±0.1 


1038.9 (fixed) 


973.7 (fixed) 


978.7 




+410±15 


+435±30 


+435 (fixed) 


+410 (fixed) 


480+50 




1135±30 


740±45 


740 (fixed) 


1135 (fixed) 


490+270 


2 S' 1 ] 


37.2±1.9 


6.0±0.4 


3.0±0.3 


10.0+0.8 


1.2+0.3 




1026.8±0.1 


1033. 1±0. 2 


1038.8 (fixed) 








+320±30 


+345±60 


+345 (fixed) 








1365±75 


935±125 


935 (fixed) 






2 s- 1 ! 


2.4±0.1 


0.5±0.07 


0.3+0.06 







SNR 0509 67.5 

P214 (LWRS) 

Acent [A] 

Vccnt [km S _1 ] 
Vfwhm [km s _1 ] 
Flux [10~ 14 ergs cm" 



1024.8+0.24 
-263+125 
3710+400 
7.8+0.6 



SNR 0519 69.0 

P214 (LWRS) 

Acent [A] 

Vccnt [km S" 1 ] 
Vfwhm [km s _1 ] 
Flux [10~ 14 ergs cm" 



2 s- 1 ! 



1026.8+0.24 
+320+70 
3130+155 
22.0+0.5 



1033.6+5.6 
+495+1630 
4975+1830 
5.4+0.36 



1039.3 (fixed) 
+495 (fixed) 
4975 (fixed) 
2.7 (fixed) 



973.5 (fixed) 
+320 (fixed) 
3130 (fixed) 
6.4+0.9 



a Rcsults quoted are for single Gaussian profile fits to each line 

b The velocity width and velocity centroid of the A1038 line are tied to the width and centroid of A1032 during 
spectral fitting. 

c The velocity width and velocity centroid of Ly 7 are fixed to those of Ly (3 
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Table 3 

Shock Parameter Estimates for LMC Balmer-Dominated SNRs 



Object 


V s [km s- 1 ] 


Plane Parallel 

(Te/T p ) 


(To/T P )o 


V s [km s- 1 ] 


Spherical Shell 

(Te/T p ) 


(To/T p )o 


SNR 0505-67.9 (DEM L 71) 

Blucshifted Shell (LWRS) 
Redshifted Shell (LWRS) 
Southern Edge (MDRS) 


775-1005 b 


<0.05 b 


1-9 


650 
900 


0.5 
0.05 


6.5 
16 


SNR 0509-67.5 

Global Spectrum (LWRS) 


5200-6300 




4-16 








SNR 0519-69.0 

Global Spectrum (LWRS) 


3200-3900 




16 


2600 


<0.2 





a When (T e /T p )o is unconstrained for plane parallel models, the shock velocities are quoted between the limits of minimal and full 
electron-proton equilibration 

b Estimatc from combined Ha and X-ray spectroscopy of Rakowski, Ghavamian & Hughes (2003) 
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DEML71 0509-67.5 




Fig. 1. — Continuum-subtracted narrow band Ha images of the four Balmer-dominated SNRs (CTIO 4m images courtesy R. C. Smith). The 
positions of the FUSE apertures (LWRS: 30" X 30" squares; MDRS: 4" X 20" rectangle) are marked. The wavelength dispersion direction 
for each of the LWRS spectra shown in Figure 2 are marked in the upper left corner of each panel. The dispersion direction for the MDRS 
pointing in DEM L 71 is perpendicular to the long dimension of the slit. 
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Fig. 2. — Combined raw FUSE spectra of the four Balmer-dominatcd SNRs. The 1A channel data (1013 A — 1067 A) are shown in the 
figure. Positions of major emission line features at the LMC systemic velocity (+275 kms - 1 ) are marked. The broad Ly /3 and O VI emission 
lines are clearly seen in the spectra of DEM L 71, SNR 0509—67.5 and SNR 0519— 69.0. A particularly strong C II absorption feature can 
be seen in the spectrum of SNR 0519—69.0 near 1037 A, arising from LMC halo absorption. There are two sets of spectra imaged on each 
detector (marked LiF 1A and SiC 1A). Each set consists of spectra from the medium resolution (MDRS, at top), high resolution (HIRS, in 
middle) and low resolution (LWRS, bottom) apertures. The object spectra are acquired through the LWRS aperture, with the remaining 
apertures containing only residual airglow. 
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Fig. 3. — Extracted LWRS spectrum of the face-on shock position in DEM L 71. The best fit line profiles are marked, along with the 
positions of major interstellar features from LMC (top tick marks) and Galactic (bottom tick marks) absorption. The model components 
included in the fit are shown in Table 2. 
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Fig. 4. — Extracted LWRS spectrum of SNR 0509—67.5. The best fit line profiles are marked, along with the positions of major interstellar 
features from LMC (top) and Galactic (bottom) absorption. The ISM absorption features have been excluded from the fit. The model 
components included in the fit are shown in Table 2. 
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Fig. 5. — Extracted LWRS spectrum of SNR 0519—69.0. The best fit line profiles are marked, along with the positions of major interstellar 
features from LMC (top) and Galactic (bottom) absorption. The ISM absorption features have been excluded from the fit. The model 
components included in the fit are shown in Table 2. 
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Fig. 6.— Extracted LWRS spectrum of SNR0548-70.4. No FUV emission is detected down to the S/N = 3 sensitivity limit of FUSE 
(~10~ 15 ergs cm -2 s -1 A -1 , at 1032 A for night only exposures. 
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Fig. 7. — The MDRS spectrum of DEM L 71, acquired from the southern rim of the SNR. The fitted emission line components (broad Ly 
(3, O VI AA1032, 1038 and Galactic + LMC H I absorption are marked. The Ly /3 airglow line (marked by the ©) is excluded from the fit. 
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Fig. 8. — The predicted broad F(Ly /3)/Fo vi(1032) flux ratios, shown as a function of the assumed electron-proton temperature equili- 
bration for the three FUV spectra dominated by edge-on blast wave emission. The measured ranges of the F(Ly /J)/Fo y/(1032) ratios are 
marked by the shaded regions. The solid and dotted lines show the predicted flux ratio curves for the maximum and minimum preshock 
neutral fractions allowing agreement between the observed and predicted flux ratios. In the case of SNR 0509—67.5 the maximum preshock 
neutral fraction (80%) is the largest plausible value given ionization conditions in the warm neutral ISM. 
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Fig. 9. — The observed and modeled Ly f} profile from the face-on spectrum of DEM L 71 (Figure 3). The solid line shows our best 
approximation to the profile from our spherical shell models (see text for details). The dashed line shows our best approximation when the 
LWRS aperture (Figure 1) extends only halfway to the edge of the shell (i.e., no edgc-on component is included). The prominent absorption 
feature and Ly f3 airglow line between —700 kms -1 and +500 kms -1 are excluded from the profile modeling. 



